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Abstract 

We consider a brane-world model in which an anisotropic brane is embedded in a 
dilatonic background. We find the background solutions and study the behaviour of 
the perturbations when the Universe evolves from an inflationary Kasner phase to a 
Minkowski phase. We calculate the massless mode spectrum, and find that it does 
not differ from that expected in standard four-dimensional cosmological models. We 
then evaluate the spectrum of both light (ultrarelativistic) and heavy (nonrelativistic) 
massive modes, and find that, at high energies, there can be a strong enhancement of the 
Kaluza-Klein spectral amplitude, which can become dominant in the total spectrum. 
The presence of the dilaton, on the contrary, decreases the relative importance of the 
massive modes. 



1 Introduction 



The brane-world scenario, developed starting from the fundamental work by Randall and 
Sundrum [1], has received enormous attention in the past years, mostly because it is pos- 
sible to generalize it to have interesting cosmological models (see [2] for a review). The 
fundamental question hence is how these models can be tested, and if there are signals 
that could allow us to distinguish between a brane-world model and a more conventional 
four-dimensional cosmological one. Studying the production of cosmological perturbations 
during inflation is a powerful tool in trying to answer to this question. Indeed much effort 
was devoted in building models (see, for example [3-11]) so as to solve the perturbation 
dynamics, which is in general very hard to deal with. 

In this work we develop a gravi-dilaton brane-world with a non-trivial dynamics on 
the brane (a lot of papers have been written generalizing the standard RS solution to a 
more general framework in which a scalar field is included, see for example [12-17]. In 
particular, following our previous work [9], we consider a higher-dimensional p-brane (with 
p > 3) coupled with a bulk dilaton. If we do not put matter on the brane, it is possible 
to solve exactly the Einstein equations to obtain a Kasner solution, with d expanding and 
n contracting dimensions, on the brane itself. The bulk equation decouples, and is solved 
by a warp AdS-like factor, while the dilaton grows logarithmically as it goes away from the 
brane. 

Turning to study the dynamics of the tensor perturbations, we find that massless and 
massive modes can be treated independently, and do not mix. We outline a procedure 
that allows us to study the production of the gravitational and of the Kaluza-Klein fluc- 
tuations; we also calculate the spectral distribution of these fluctuations, amplified during 
an imaginary phase transition from the inflationary Kasner regime to a simple Minkowski 
phase. This is done for the massless mode, and for both light and heavy-massive modes 
with respect to the curvature scale at the transition epoch. Our results show, differently 
from what happens in other models present in the literature, that if the transition occurs 
at high enough energy, there can be a strong enhancement of the contribution of the KK 
modes; this should be in principle observable, since, in this regime, the complete spectrum is 
found to be quite different from what would be expected if there were only the conventional 
4D graviton. We also find that, on the contrary, the dilaton effect results in lowering the 
amplitude of the KK perturbations. 

The paper is organized as follows: in section [21 we present and solve the background 
equations. In section |3] we study the perturbation around this background, and define the 
action that controls the dynamics of each single mode (massless and massive) on the brane. 
Section 0] outlines how to develop the canonical analysis of the mode action, and how to 
obtain the correct canonical equation that describes the evolution of the mode during the 
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phase transition. In section |S] we obtain the correct expression for the spectral amplitude, 
which is specialized for the graviton in section|BJ and for the KK modes in section[7| Finally, 
in section|Sl we make some considerations on the results obtained, and draw our conclusions. 



2 Background solutions 

Let us consider a model in which a brane is non-minimally coupled to a bulk dilaton eft. 
We work in a -D-dimensional space (D = p + 2), and set the brane fixed at the origin: 
X D ~ l = z = 0. The action is 
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(2.1) 



where jab = 9AB + tiatib is the induced metric on the brane, T p is the brane tension, and 
L£> is the bulk cosmological constant. 

Variation of this action with respect to the metric and to the dilaton gives the Einstein 
equations: 
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and the dilaton equation: 
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These equations, written as above, already include the Israel junction conditions [18], which 
can be deduced by integrating along a small interval across the brane on a geodesic perpen- 
dicular to the brane itself [2, 19] (see also [20]). They read, for the Einstein equations (|2.2j) 



[k 



(2.4) 



where the square brackets denote the difference between the left and the right limiting value 
on the brane: 

[/] = ( lim - lim ) /(e), (2.5) 



lim — lim I f(e), 

and Jzab is the extrinsic curvature on the brane. 

To solve the equations we set the following ansatz on the metric: 



ds 2 = f 2 {z) (dt 2 - a 2 (t)dx 2 - 6 2 (i)dy 2 - dz 2 ) , 



(2.6) 
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i.e. we allow the metric to be anisotropic, and impose that d spatial dimensions expand 
and the other n contract (p = d + n). We also impose that the dilaton depend only on the 
extra-dimension: <j) = <p(z). The equations specialize in: 
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(2.10) 
(2.11) 



where dots (primes) denote differentiation w.r.t. t (z), and H = a/a, G = b/b, F = f'/f. 
Note that the singular part is absent from the above equations, because it will be taken into 
account by satisfying the Israel junction conditions. Inserting the ansatz Q2.6J1 into (|2.4|) we 
get the metric junction condition: 
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The dilaton junction condition is obtained directly by applying on eq. (|2,3j) the integration 
procedure described above, as in standard one-dimensional quantum mechanics, and it reads 
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The Einstein equations (|2.7H2.1()|) can be decoupled, and the time dependence is solved, 
as in [9], by the Kasner solution: 



a(t) 
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(2.14) 
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which can describe a superinflationary solution [21,22] on the negative branch of the time 
axis, if we choose the minus sign in I and allow for 2 or more internal dimensions (of course, 
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in the framework of superstring theory there is room for up to 5 internal compact dimensions 
if we want 3 external and one warped). So we are left with the z dependent part of the 
Einstein equations, which can be rearranged as: 

- P F'-p 2 F 2 = 2L D e ai(f> f 2 
-pF'+pF 2 = <j)' 2 , (2.15) 

and the dilaton equation (|2.11j) . which nevertheless depends on the other equations by 
means of the Bianchi identities, as expected. 

To solve eqs. (|2.15(l we seek for a solution of the form [14,23]: 

/(*) = i i +£ a 

<j>(z) = logfl + iQ ; (2.16) 

here z$ is a positive constant (which corresponds to the AdS length in the usual RS model), 
and the solutions are intended to be in the z > region, since we consider a Z2 symmetric 
background. By inserting these expressions in (|2.15j) and taking into account the Israel 
junction conditions ([2.12ft and (|2.13ft we get the following expressions for the exponents 

paf — 4 paf — 4 

and the following relations between the parameters 

2(p + l) 
T p 8p 
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To get a positive brane tension and a negative bulk cosmological constant, we assume that: 

-4= < cn < 4- (2-19) 

From (|2.18jl it is easy to get the relation between the brane tension and the bulk cosmological 
constant 

M~p~ = 4 ^/ 4 £±i-\ 2 ' (2 " 20) 
which is the generalization of the fine-tuning relation in the standard RS scenario. 
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It is possible to obtain a different class of solutions, by imposing a different ansatz: 



/(*) = e 



-az/z 



m = -p-. (2.21) 

This solution saturates the bound (|2.19j) a± = 2/y/p, and has the exponents related as 
follows 

(3 = ±y/pa. (2.22) 
The relation between the tension and the cosmological constant is unchanged 



w = (2 - 23) 

Nevertheless, we will only consider in what follows the first kind of solution described above. 

3 Perturbation equations 

In this section we are going to derive the equations of motion for the tensor perturbations of 
the metric, S^gAB = h>AB- We set the dilaton perturbation equal to zero, 5^<j) = 0, because 
it would decouple from the tensor fluctuations, impose that the perturbation depends on 
only the external spatial dimensions and work in the transverse-traceless gauge: 

hoA = Ka = h D -i,A = 0, g lj hij = Vjh { J = 0. (3.1) 

The second-order perturbation of the countervariant indices of the metric and of the dilaton 
is set to zero as well. To obtain the expression of the induced metric, we use the equation 

gAB n A n B = -1 = (gf B + h AB ) (n^ A + ^V) (n^ B + S^n B ) . (3.2) 

The zeroth-order part of this expression is equal to —1, and we can ignore orders higher 
than the first. Then eq. (|3,2j) becomes: 

2g AB n^ A 5Wn B = = -|^n 5 . (3.3) 

So the extra-dimensional part of the normal unit vector vanishes, and therefore the complete 
vector is left unchanged as well, 5^n A = 0. This shows that the induced perturbed metric 
is: 

5 (1) 7ab = hAB = h A B\ z=0 ■ (3.4) 

The perturbation equations can be obtained by variating the action (|2.1j) perturbed to 
order h 2 . After some algebraic manipulation and making use of the background equations 
(I2.7H2.10B . we find the final form: 

M d 



5^S {a) = — / d d+1 xdz a d b n f p 



■ 2 (dkh (a) ) 2 



k=l 



(3.5) 
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where we have assumed that the internal dimensions have been compactified on a com- 
pact manifold of size M~ n . This is the action for the single polarization mode: S&'S = 
Yl( a ) <^ 2 ^(a)- The polarization mode is defined via the spin-2 polarization tensors, 

25 ab . From now on we will omit the 



hi 



h(a) e ij\ which satisfy the relation: e'^ 'e { 



l (a) k (b) 



polarization index a. Variation of (|3.5|) leads to the equation of motion for each mode of 
the tensor perturbations, which is, as expected, the D-dimensional covariant d'Alembert 
operator on the background considered: 

V 2 



h + (dH + nG) h T h - h" - pFti = 0. 



(3.6) 



1 T P c (p+l)a2(, 



We now need a prescription to project the perturbation equations, which is free of 
singularities, on the brane. This condition can be obtained by perturbing to the first order 
the Israel junction condition (|2.4|) . Making use of the particular form of the metric Q2.6JI 
we derive the condition: 

z =o = (/ 

To solve the perturbation equation (|3.6j) . we can expand its solutions on an orthogonal 
basis of functions of z (the precise form of the scalar product, and consequently of the 
normalization of the functions, will be discussed shortly): 

h(t, x l , z) = v (t, x l )ip (z) + dm v m (t, x % )ij) m (z) = h + dm h 
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(ti + 2Fh) 
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(3.8) 

Inserting this expansion in (|3.6|) and in 1)3 ,7|) leads respectively to the two equations (which 
are valid for both the massless and the massive modes): 



v m + (dH + nG) v, 
and to the condition 

dm v m 
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(3.9) 
(3.10) 
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Let us now substitute the decomposition (|3.8j) into the action (|3.5|) . First we consider 
the zero mode. After integrating by parts and making use of 1)3. 9 J) we get 

(^ ) 2 " 



M° 
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So it is clear that the auxiliary field tpQ 
be normalized to unity: 



Mf 2 -00 has the correct canonical dimension to 



dz \^ Q {z)\ 2 = M / dz f p \Mz)\ 



(3.13) 
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From this we can write down the final form for the action that describes the evolution of 
the zero mode 

d 



8® S 



M d ~ 
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4|Vo(0)| s 
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and the effective Planck mass of the 4-dimensional effective graviton can be easily read off 
to be: 



d-l 



M d ~ 



(3.15) 



4|Vo(0)| 2 ' 
as in the ordinary RS scenario [1,9]. 

Now we turn to study the massive modes. Again, after a little algebraic manipulation, 
we arrive at 

rd 
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In this case the auxiliary fields with the correct canonical dimension are i/j m = f 2 ip m , which 
satisfy a Schrodinger-like equation 

—F 2 + P -i 
4 2 

therefore, they can be normalized as a complete set of orthonormal solutions 



m 



$m = 0; 



(3.17) 



dz i> m 4> m > = dz f p ^ m ^ m ' =5{m-m!). 



(3.18) 



By making use of this relation in (|3.16|) we finally find that the action (|3.5[) can be thought 
of as the sum of the massless mode action with an infinite set of massive actions, one for 
each single massive mode: 



S (2) S = 6 (2) Sq + J dm 6 (2) Sm _ 



(3.19) 



The massless mode action 8^ So of (|3.14|) can be interpreted as the 4-dimensional graviton, 
while the form of the massive mode action is: 
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From (|3.2U|) we can deduce the value of the generalized "Planck mass" for each mode as 
the coefficient that multiplies the mode action: 

M ™ = rar (3 - 21) 

At this stage, two remarks are in order. The first is that the massive mode action 
does not contain terms that couple the massive modes to each other, i.e. each KK mode 
is "free" with respect to the interaction with the others. This is due to the fact that the 
perturbation equation Q3.6JI decouples (at least in the first-order approximation), so that we 
can treat the two equations (|3.1U|) independently, and the Israel junction condition ()3.11j) 
can be satisfied for each massive mode independently, as we will do in section (this is not 
true in general, see for example [10] in which a similar method leads to an action for the 
modes that contains coupling terms). Moreover, it is worth stressing that the massive mode 
action is not adimensional, since it has the dimension of a length, [6^S m ] = [M -1 ]. This 
is of course a consequence of the fact that the spectrum of the KK modes is continuous, 
and indicates that physical quantities will be obtained by integrating over a suitable mass 
interval. But, formally, for what concerns the manipulation to obtain the massive mode 
spectrum, we will treat ()3.2U|) as a genuine canonical action. Nevertheless, the important 
result here is that modes do not mix and can be treated separately. 

4 Canonical analysis of the perturbed action 

To calculate the spectrum of the fluctuation, we need to put the actions (j3.14|) and (j3.2Uf) 
in a canonical form. In order to do this we will adopt the conformal time to describe the 
evolution of the system (from now on a dot will denote a derivation with respect to rj): 

dv = ^r v (4.1) 
a{t) 

and introduce the pump field 



M d ~ l 1 d-ln , , M d 1 d-ln y . 

^ ) = V^v^f^< u.) = V— vwoT™ (42) 

Through the pump field it is possible to introduce the canonical field 1 : 

u m (rj, x 1 ) = € m (rj)h m (r], x l ). (4.3) 

1 We stress again that the canonical field for the massive modes actually does not have a canonical 
dimension. 
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In terms of this canonical field, the action can be expressed in a normal form (valid both 
for m = and for m/0): 



S^Sm = 5 I d".r<lil 



U 2 



Era \2 i | Cm 2 2| 2 

(Oi«m) + 7 ma \u m 

i=l \ u / 



(4.4) 



This form for the action makes it possible to adopt the standard quantization procedure: 
we promote u m to operators and impose the canonical commutation relations among the 
fields and their conjugate momenta 7r m (r],x l ) = u m (rj,x l ): 

[7x0(77, x), 7^(77, x')] = [710(77, x), 710(77, x')] =0 
[u m (r],x),u m/ (rj,x f )] = [71^(77, x), 7r m / (77, x')] =0 
[7x0(77, x), 710(77, x')] =i(5 d (x'-x) 

[7x^77, x),7r m / (77, x')] = i5(m - m')5 d (x - x.'); (4.5) 

we then express them in terms of their Fourier components: 

d d k 



u m {r),x l 



/d k r 
(2n)d/2 l u m,k(v)a m (k)e l]l -* + <, >fc (f7)<4( fc ) e ~ 



tkx 



d d k 



-tkx 



(2n) d / 2 

The operators {a m (k)} can be made to obey a canonical oscillator algebra 

[ao(k),a (k')] = [a (k),a (k')] =0 
[a m (k),a m ,(k')] = [a^(k),a^(kO] =0 
a (k),4(k')] = <5 d (k-k') 
a m (k),a^(k')l =5(m-m')^(k-k'), 



(4.6) 



(4.7) 



and therefore they can be interpreted as a set of creation-annihilation operators, provided 
the wavefunctions {u m ^{f])} satisfy the equation (which has, again, a Schrodinger-like form): 



where we have denoted 



v m ,k(v) = \ k 2 + m 2 a 2 (r]) - 



&n (??) 
im{v) 



k 2 - V mtk (r]), 



and they are normalized according to the wronskian condition (as expected) 

u m ,k{rj)u* m!k {ri) - 74^(77)7^(77) = i. 



(4.8) 



(4.9) 



(4.10) 
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We should stress that, even if the massive fields u m (with m / 0) do not have the 
correct canonical dimension, [u m ] = [M^~], this discrepancy is compensated by the di- 
mension of the creation-annihilation operators, as can be deduced from (|4.7|) . So the mode 
coefficients u m ^ still have the proper canonical dimension of quantum-mechanical wave- 
functions: [u m> jt] = [M - 5]. This will be important when we will need to normalize the 
initial fluctuations, as we will do in the next section. 

5 Derivation of the spectrum 

The aim of this section is to calculate the spectral distribution for both the massless and 
the massive gravitational modes, using the actions developed in the last section. Hence, let 
us imagine a phase transition in which the Universe evolves from the inflationary regime 
described by the Kasner solution presented in section [2] to a simple final Minkowski era. 
Since we are interested in models that have 3 external dimensions, from now on we will 
assume d = 3. Moreover, to have an inflationary expanding initial phase, we need, as 
observed before, 2 or more internal dimensions. For the sake of simplicity we will take 
the value n = 2, but the consideration we will make can easily be extended to a different 
number of internal dimensions. The phase transition will be analysed by making use of 
the sudden transition approximation (see, for example, [24]), so that the geometry changes 
instantaneously at the conformal time —rjx. This approximation gives reliable results only 
if the frequency of the amplified modes is much lower than the transition velocity, which 
is estimated by the curvature scale at the end of inflation H\ ~ 1/niai = fci/oi = k% (of 
course a± = a(— r/i) = 1). For this reason, Hi represent a cut-off frequency. 

In this framework, the Bogoliubov coefficients that describe the transformation from 
the | in) to the \out) states can be obtained by simply imposing continuity of the Fourier 
coefficients u m ^ and their time derivative at the transition time This procedure is 

nevertheless effective only if we have an unambiguous way to normalize the \in) functions 
to pure positive norm states. This is indeed what happens in our case, since it is easy to 
see that the background solutions, written in conformal time, are 




with the parameters I and [i given in (|2.14|) : since, if we want an inflationary solution, it 
must be I < 0, we see that a(rj) vanish as rj — ► — oo. On the other side, the expression for 
the pump field during the Kasner regime is (note that it is independent from the number 
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of internal and external dimensions): 



^-s—mri- ui " ) = i-^mtb (5 - 2) 

but of course the ratio £ m /£; m always has the same behaviour, like r\~ 2 . So the whole 
potential V m ^(r]) goes to zero at the infinite past, and in this limit the mode equation 
reduces to that of a massless non-interacting field, Um,k+k 2 u mj k = 0. Moreover, as explained 
in section |1J the dimension of the mode coefficient is the same as expected, so we can guess 
the correct initial expression for it: 



u m , k ( V ) = ^=e~^. (5.3) 



1 

2k 

With this in hand, what we need to do is to find a complete solution for ()4.8j) fixing the 
free parameters in a suitable way to match the asymptotic solution ()5.3j) . and then to solve 
the linear system 

U ml(.-Vi) = a m (k)u+ i ° ut \-i 11 ) + p m (k)u^ { ° ut) (-m) 



u 



S(-r/i) = c^{k)tiffi\- m ) + PUk)n^ Ut \-Vi), (5-4) 

to find the Bogoliubov coefficient f3 m (k). Here u^°^ and ujfu' are respectively the 
positive- and negative-frequency eigenfunctions in the Minkowski phase. The Bogoliubov 
coefficient represents the amount of particles created by the gravitational fields in the mass 
interval [m,m + dm]. From this we can obtain the energy density per logarithmic interval 
unit for each mode: 

di^k = T* mk) \ (5 - 5) 

= ^Vk^+^\P m (k)\ 2 dm. (5.6) 
a log ft vr 

These expressions should be normalized with their respective coupling constants Q3.15|) and 
(|3.21[) . Note that the different dimensions of the spectral distributions are balanced by 
the different dimensions of the massless and massive coupling constants, which explains the 
quantization procedure with the "odd" dimensions of section 0J Then the massive spectrum 
should be integrated over all masses (eventually one can be interested in the contribution of 
a particular mass interval). To obtain a dimensionless quantity, we then choose to normalize 
the spectral distributions just obtained to the scale curvature H\ at the end of inflation. 
The spectral distribution can therefore be cast in its final form: 

n(k) = n (k) + n KK (k), (5.7) 



n 



with 



Qo(k) 



k 4 



tt 2 H 2 M 2 1 
Okk(^) = / dm £l m (k) 



dm 



-Vk* + m2\f3 m (k)\ 2 . 



(5.8) 
(5.9) 



In the next sections we will apply this procedure to derive explicitly the spectral amplitudes. 



6 Spectral distribution for the massless mode 

The equation for the Fourier modes (|4.8j) becomes, in the massless case: 

1 



uo,k + 



u ,k + k 2 u 0:k = 7] > -771, (6.1) 

so the I in), solution which asymptotically tends to (|5.3|) . and the \out) solutions are: 



(in) 
0.A; 



7T 



j\v\H^(k V ) 

-ikr) 



II: 



+(out) 
0,k 



II, 



-(out) 
0,k 



2k 

Jkr] 



2k 



(6.2) 



(2) 

where Hq is the Hankel function of the second kind. This gives, for the Bogoliubov 
coefficient (3, the expression 



A) (*0 



7T 



Jkrji 



2ki 11 H^\ki 11 ) + i ( 2H ( Z> (k m )km - H ( z> {k Vl 



(2)', 



(6.3) 



32/cni 

The expression in square brackets can now be approximated for small values of the argument 
of the Hankel functions, since our approximation is valid if kr/i <C 1, and the leading term 
is logarithmic. So the energy density (J5.5|) takes the form 2 



d PQ _ H l fk\ log2 k_ 
d\og k 87r 3 \k\) k\ 



(6.4) 



The next step is to solve eq. (j3.9j) for the massless mode. Its general solution can be 
written as 3 



ijj (z) = c + Ci 



1 — oa \ zq 



l-5a 



(6.5) 



2 In plotting the spectrum we do not need to use the approximation. Actually, in the forthcoming 
calculations, we heavily rely on numerical computations. 

3 We recall that, in what follows, the numerical values are only valid in d = 3 and n — 2 dimensions. 
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but the normalization condition (|3,13|) imposes c\ = 0, so the only acceptable solution is 
the constant one tpo(z) = cq. The free parameter cq can be calculated using (|3.13|) as well, 
so we get: 

Finally, the expression for (|5.8|) is 



1_ ( Hi\ 2 ( k_\ \ > /• 



Mr 



87T 3 V M P / V^l 

z M 3 



p 2(-l-5a)' 

and it is plotted in Fig. ^for different values of the parameters ZqHi and a. 
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Figure 1: The spectral amplitude of the massless mode. In a) the spectrum is evaluated at 
a = — 1 and for zqH\ = 7 (dashed line), ZqHi = 20 (dotted line) and zqHi = 50 (dot-dashed 
line). In b) the spectrum is evaluated at zqHi = —7 and for a = — 1 (dashed line), a = —5 
(dotted line) and a = — 8 (dot-dashed line). 



7 The KK spectrum 

First of all we will evaluate the massive coupling constant. In order to do this we have to 
solve eq. (|3.17|) for ra/0. Using the background solutions (|2.16l) . it is easy to see that this 
equation is again a Bessel-like equation of the form: 

2 5a(5a - 2) 1 



< + ( m 2 - — 4 >- ^-—y )^ m = 0, (7.1.) 
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and its general solution can be written as 

ipm(z) = y/z + ZQ [cx(m)J v (m{z + z )) + c 2 {m)Y v (m(z + z ))] , (7.2) 

where J v and Y v are of course the Bessel functions of the first and of the second kind, and 
the parameter v is related to the dilaton coupling constant by the relation: 

1 — 5a 



, (T.3) 

The constants c\ and C2 are fixed by using the junction condition ([.3.11)1 and the normaliza- 
tion condition (|3.18|) . Prom the junction condition we get, using some algebraic properties 
of Bessel functions 

J u -i(mz ) 

c 2 = -tp rci. (7.4) 

From the second one, using their orthonormality relation 4 , we get 

ci(m) = ym (7.5) 

y Ju-i( mz o) + Y u-i( mz o) 

Finally we can write the properly normalized solution of ([7.1)1 as 

? <\ r~i — ! ?[Y 1/ -i(mz )Jv(m(z + zo))- Jv-i(mz )Y u (m(z + zo))] 

Wm{z) = ^Jm(z + zq) = . (7.6) 

^Jl_ x {mzQ) +Y^ 1 (mz ) 

Now we turn to consider eq. ([4.8)1 . Actually it is very hard to solve it in its complete 
form, so the best we can do is to seek for particular limits in which the equation simplifies 
a little. As already pointed out, in the phase transition we considered, an energy scale 
naturally emerges, i.e. the curvature scale at the end of the inflationary epoch. In order to 
have a complete understanding of the effects of the KK modes, we will study what happens 
to lighter or heavier modes as regards to this curvature scale. Firstly, we limit our attention 
to modes that are lighter than the curvature scale at the end of inflation, m <C H\. In 
this case we can neglect the mass term in ()4.9(l . so that the mode equation reduces to the 
same as obtained for the massless mode, eq. (|6.1j) . It is easy to see that even if the \out) 
solutions are massive waves, the spectral distribution remains unchanged, as in eq. (|6.4|) . 
So, the only contribution to the integral, which must be calculated from m = to m = H\, 
comes from the massive coupling constant. The integration leads to a change of the mass 
parameter that controls the normalization of the spectral amplitude, while the shape of the 
spectrum remains unchanged. We obtain 

niightC*) = dm n m (k) = ± (^-Y ( £Y log 2 A, (7.7) 



, , , ^ 2 

2vr 3 iMj \hj h' 



/o Ml,) 4 \J 



M t = ( / W ) = — ( I dm |iM0)|' ) • (7.8) 



See for example the classical electrodynamics textbook by J.D. Jackson 
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In Fig. [21 we show a numerical estimate of the behaviour of the spectrum for different values 
of the parameters. 
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Figure 2: The spectral amplitude of light KK modes. As before, in a) the spectrum is 
evaluated at a = — 1 and for zqH\ = 7 (dashed line), ZqH\ = 20 (dotted line) and zqH\ = 50 
(dot-dashed line). The three curves overlap, so that they are quite confused. In b) the 
spectrum is evaluated at zqH\ = —7 and for a = — 1 (dashed line), a = — 5 (dotted line) 
and a = —8 (dot-dashed line). 



Next, we try to solve (|4.8[) for the heavy modes. To do this, we use a WKB-like 
approximation for the mode function u m ^ [25,26], which can be written as: 



(in) 



exp 



i / drfu m>k {rf) 



(7.9) 



This approximation is valid if the variation of the frequency is small with respect to the 
frequency itself. More precisely, (see, for example, [27]), one must have 



_ 3 Um,k 
m.fe 



1 Q, 



m,k 



2w 3 , 

m.k 



< 1. 



(7.10) 



Since the parameter e grows monotonically in time, we only need to check where our ap- 
proximation is valid at the transition epoch rji, and we are guaranteed that it always holds 
before that time 5 . It is not difficult to see that |e(— ryi)| <C 1 if ma\ S> fei, which is exactly 
the regime we wish to investigate. This is not an unexpected result, since it is known [27] 



3 To be more precise, e actually reaches a maximum, but well after the transition time. 
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that this approximation works better if the transition is adiabatic, so that the variation of 
the curvature is slow with respect to the energy. But this means that the sudden transition 
approximation is no longer valid. Nevertheless, we could still get clues on the form of the 
spectrum by analytically continuing the function LO m ^ and evaluating the integral in (|7.9f) 
on a suitable path. We consider a semicircle on the upper half-plane. Its radius should be 
chosen to be greater than 771 , to stay in the region in which the WKB approximation is 
valid, but not too big, since we expect the mass term to be dominant with respect to the 
pure frequency term k 2 . In this regime the function uj m ^ can be approximated as 

k 2 1 £ m 

w m , fc (r?) ~ 7710(77) + — - r-r-f 1 , (7.11) 

2ma[rj) 2ma(rj) £ m 

and has, as one could expect, the form of non-relativistic energy for a massive particle, 
plus a term that describes the direct interaction between the particle and the background 
geometry, which however is negligible. Note that, in the approximation (|7.11|) . uj m! k has 
only a singularity at the origin, so we can actually shrink the radius of the integration path 
till R ~ 771 without changing the result of the integration [26] . It is not difficult to see that 
the square modulus of the Bogoliubov coefficient is given by: 



\P m (k)\ 2 = exp 



2 Im (J dr] uj m>k (r]) 



(7.12) 



The integral is easy to solve if we make use of the approximation (|7.11|) . so we get: 



\t3 m {k)\ 2 = exp 



mai k k 



2sin7rA - — ; =- — ; =- + 



fei(l + A) 2moiA;i(l-A) 8mai(l + A) 



m > H x . 
(7.13) 

This expression can now be substituted in (|5.9|) to obtain the final form for the Kaluza-Klein 
spectrum. The integral in this case runs from m = H\ to infinity. 



-DC 



A: 3 



fibeavy(*0 = / dm 2ir2A/fS V k 2 + m 2 \(5 m (fc) 1 2 ■ (7.14) 

Of course this integration can only be carried out numerically. The result is presented in 
Fig. |31 for the same values of the parameters zqH\ and a used in the previous computations. 



8 Comments and conclusions 

In this paper we have presented a simple model in which a brane coupled with a bulk 
dilaton evolves from an initial inflationary Kasner phase to a final Minkowski era. Then 
we have studied the evolution of the tensor perturbations on this background, paying par- 
ticular attention to their normalization to an initial state of vacuum fluctuation. We have 
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Figure 3: The spectral amplitude of the heavy KK modes. As in the other figures, plot a) 
represents the spectrum evaluated at a = — 1 and for ZoHi = 7 (dashed line), zoHi = 20 
(dotted line) and zqH\ = 50 (dot-dashed line). Again, overlapping of the curves do not 
permit to distinguish them. Plot b) represents the spectrum evaluated at zqH\ = —7 and 
for a = — 1 (dashed line), a = — 5 (dotted line) and a = —8 (dot-dashed line) 



found that the massless and the massive modes can be treated independently, and we have 

evaluated (analytically where it was possible, numerically elsewhere) the spectrum of the 

tensor perturbations for the massless mode (to be identified with the graviton) and for both 

the ultrarelativistic and the non-relativistic massive modes. Of course the total spectrum 

an observer would detect should be approximatively the sum of the three. 

The behaviour of the plot we presented suggests some comments: As widely expected, 

the relative importance of the KK corrections on the total spectrum grows as the AdS 

curvature increase 6 , because [4,5,19] deviations from classical general relativity become 

more and more relevant as the energy increase. On the other hand, the coupling between 

the dilaton and the brane has the opposite effect. In fact, as u\ approaches its limiting 

value 2/y/p (see eq. (|2.19)l ). the massive perturbations are highly suppressed. We can have 

a better understanding of the relative importance of the contribution of the KK modes by 

6 Note that we have chosen to normalize the spectrum to the end of inflation curvature, so that an increase 
of the AdS curvature actually corresponds to an increase of the curvature scale at which the transition occurs. 
This is the reason why it seems that the massless mode spectrum is influenced by the change in zq and the 
KK spectrum is not. 
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defining the ratio: 

r _ ^KK ^ ^ 

£l + ^KK ' 

and by plotting its behaviour with respect to the two parameters we are interested in. This 
is done in Fig. ^ 
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Figure 4: Plot a) shows the behaviour of r with respect to ZoH% for the three values of 
a used elsewhere: a = — 1 (dashed line), a = — 5 (dotted line) and a = —8 (dot-dashed 
line). Plot b) is complementary to the first one, and shows the behaviour of r with respect 
to a for the three different values of zqHi = 7 (dashed line), zqHi = 20 (dotted line) and 
zqHi = 50 (dot-dashed line). 



As we can see, the KK contribution is very low at small values of zqH\] however, it 
rapidly increases, to become completely dominant in the high energy regime, differently 
from what happens in other models present in the literature [5-7, 11]. Of course, this 
effect could in principle be relevant in observational tests concerning the amplification of 
tensor perturbation during inflation. Actually, since the spectrum obtained in our simple 
evolution model increases with the frequency, it is not difficult to tune the free parameters 
so as to satisfy COBE and pulsar constraints. Moreover, the relevance of the massive 
modes contribution can be lowered by the presence of the dilaton. In fact, at relatively high 
energies, the effect of the bulk gravitons can be very strong if the dilaton is absent or weakly 
coupled, but it becomes negligible as the coupling increases. This indicates that, even at 
high energies, the enhancement of the amplitude of the KK spectrum can be "cured" by 
a suitable choice of coupling parameter, leading to a spectrum that is quite similar to the 
4-dimensional one. 
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The decisive test for models with strong production of KK modes will be the future 
detection of stochastic gravitational waves at gravitational antennae. For this purpose, 
more accurate and realistic models are needed, in which the presence of matter on the 
brane is taken into account, and a smooth transition from the inflationary to the standard 
cosmological evolution is considered. This will be the subject of forthcoming papers to 
investigate which features of the simple model under discussion can be generalized to these 
more realistic models. 
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